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ABSTRACT 

We present a model for the formation of massive black holes (~ 1000M Q ) due to stellar-dynamical 
processes in the first stellar clusters formed at early cosmic times (z ~ 10 — 20). These black holes are 
likely candidates as seeds for the supermassive black holes detected in quasars and nearby quiescent 
galaxies. The high redshift black hole seeds form as a result of multiple successive instabilities that 
occur in low metallicity (Z ~ 1O~ 5 Z ) protogalaxies. We focus on relatively massive halos at high 
redshift (T vir > 10 4 K, z > 10) after the very first stars in the Universe have completed their 
evolution. This set of assumptions ensures that (i) atomic hydrogen cooling can contribute to the gas 
cooling process, (ii) a UV field has been created by the first stars, and (iii) the gas inside the halo 
has been mildly polluted by the first metals. The second condition implies that at low density H 2 is 
dissociated and does not contribute to cooling. The third condition sets a minimum threshold density 
for fragmentation, so that stars form efficiently only in the very inner core of the protogalaxy. Within 
this core, very compact stellar clusters form. The typical star cluster masses are of order 10 5 M© 
and the typical half mass radii ~ 1 pc. A large fraction of these very dense clusters undergo core 
collapse before stars are able to complete stellar evolution. Runaway star-star collisions eventually 
lead to the formation of a very massive star, leaving behind a massive black hole remnant. Clusters 
unstable to runaway collisions are always the first, less massive ones that form. As the metallicity of 
the Universe increases, the critical density for fragmentation decreases and stars start to form in the 
entire protogalactic disk so that i) accretion of gas in the centre is no longer efficient and ii) the core 
collapse timescale increases. Typically a fraction ~ 0.05 of protogalaxies at z ~ 10 — 20 form black hole 
seeds, with masses ~ 1000 — 2000M©, leading to a mass density in seeds of a few ~ 1O 2 M0/Mpc -3 . 
This density allows enough room for black hole growth by accretion during the quasar epoch. 
Subject headings: black hole physics - instabilities - stellar dynamics - galaxies: nuclei - galax- 
ies:formation 



1. INTRODUCTION 

Supermassive black holes (BHs) are routinely detected 
in the centre of galaxies, both in nearby quiescent galax- 
ies, and as the engines that power quasars and active 
galactic nuclei. Observationally, we can trace quasars 
until very high redshifts. Luminous quasars are detected 
in the Sloan survey (e.g. Fan et al. 2001) at z ~ 6, 
corresponding to a time when the Universe was not even 
one billion years old. This implies that the first BH seeds 
must have formed earlier on. The currently favored sce- 
nario for BH seed formation relies on the remnants of the 
very first generation of metal-free stars (Population III, 
PopIII) . The first stars are believed to form at z ~ 20 in 
halos which represent rare high peaks of the primordial 
density field. Simulations of the fragmentation of zero- 
metallicity protogalaxies suggest a very top-heavy initial 
stellar mass function, and in particular the production 
of very massive stars with mass > 100M Q . PopIII stars 
in the mass range 140 < m* < 260 M Q are predicted 
to make pair-instability supernovae. If zero metallicity 
very massive stars form above 260 M Q , they will rapidly 
collapse to BHs with little mass loss (Fryer et al. 2001), 
leaving behind BHs with masses ~ 10 2 M© (Madau & 
Rees 2001). 
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Although this path to BH seed formation seems very 
natural, large uncertainties exist on the final mass of 
PopIII stars. Even recent simulations (Gao et al. 2006) 
have not clarified if PopIII stars are indeed very mas- 
sive, and in particular if they are above the threshold 
(~ 26OM ) for BH formation (but see Freese et al. 2008, 
Spolyar et al. 2008, Natarajan et al. 2008). Further- 
more, metal abundances in extremely metal poor Galac- 
tic halo stars, which are commonly thought to trace the 
enrichment products of the first generation of stars, are 
incompatible with the yield patterns of zero metallicity 
very massive stars (> 100M Q , Tumlinson et al. 2004). 

Alternative routes to BH seed formation have been ex- 
plored (Haehnelt & Rees 1993, Loeb & Rasio 1994, Eisen- 
stein & Loeb 1995, Bromm & Loeb 2003, Koushiappas 
et al. 2004, Begelman, Volonteri & Rees 2006, Lodato 
& Natarajan 2006), typically exploiting gas-dynamical 
processes in metal-free galaxies. Gravitational instabil- 
ities can indeed lead to a vigorous gas inflow into the 
very central region, supplying the necessary matter for 
the formation of BH seeds. The typical conditions that 
lead to efficient gas infall and BH seed formation can be 
summarized as (i) the host is massive enough that the 
virial temperature exceeds T V - 1T > 10 4 K, so that gas is 
able to cool and collapse via atomic hydrogen cooling, 
(ii) Molecular hydrogen does not form as the gas cools 
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and condenses, (iii) The gas has primordial composition, 
so that metal line cooling is non-existent. 

However, massive halos (T v ; r > 10 4 K, masses ~ 
10 7 M Q ) are likely built-up from mini- halos (T v j r < 10 4 
K) that had collapsed earlier on. Some of these ha- 
los might have experienced PopIII star formation, and 
should be enriched with at least some trace amount of 
metals. Fragmentation and formation of low mass stars 
starts as soon as gas is polluted by metals created in the 
first PopIII stars. Efficient gas collapse, leading to BH 
seed formation, is mutually exclusive with star forma- 
tion, as competition for the gas supply limits the mass 
available. 

However this first episode of efficient star formation 
can foster the formation of very compact nuclear star 
clusters (Clark et al. 2008, Schneider et al. 2006) where 
star collisions can lead to the formation of a supcrmassive 
star, possibly leaving a BH remnant with mass in the 
range - I0 2 - I0 4 M Q (Omukai et al. 2008). 

The possibility that an "intermediate-mass" BH could 
in principle form as a result of dynamical interactions 
in dense stellar systems is a long standing idea (Begel- 
man and Rees 1978; Freitag et al. 2006b, a, Ebisuzaki 
et al. 2001; Portegies Zwart and McMillan 2002; Miller 
& Hamilton 2002, Gurkan et al. 2004). During their 
lifetime collisional stellar systems evolve as a result of 
dynamical interactions. In an equal mass system the 
central cluster core initially contracts as the system at- 
tempts to reach a state of thermal equilibrium: energy 
conservation leads to a decrease in the core radius as 
evaporation of the less bound stars proceeds. As a re- 
sult the central density increases and the central relax- 
ation time decreases. The core then decouples thermally 
from the outer region of the cluster. Core collapse then 
speeds up as it is driven by energy transfer from the cen- 
tral denser region (Spitzer 1987). This phenomenon is 
greatly enhanced in multi-mass systems like realistic star 
clusters. In this case the gravothermal collapse happens 
on a shorter timescale as dynamical friction causes the 
more massive stars of mass m, to segregate in the centre 
on a time-scale t^i — ((m)/m)t r h (where t r h is the half 
mass relaxation timescale and (m) is the mean stellar 
mass in the cluster). If mass segregation sets in before 
the more massive stars evolve out of the main sequence 
(<~ 3 Myr) , then a sub-system decoupled from the rest of 
the cluster can form, where star-star collisions can take 
place in a runaway fashion that ultimately lead to the 
growth of a very massive star (VMS) (Portegies Zwart et 
al. 1999). 

Yungelson et al. (2008) study the fate of solar com- 
position VMSs in the mass range 60-1000 M . They 
find that all VMSs are likely to shed most of their mass 
via winds well before experiencing a supernova explo- 
sion. Solar composition VMSs are therefore expected to 
end their lives as objects less massive than ~ 150M Q , 
collapsing into BHs with mass < 15OM or exploding as 
pair- instability supernovae. The growth of a VMS should 
be much more efficient at low metallicity. Low metallic- 
ity can modify the picture in different ways. First, at 
sub-solar (but still not primordial) metallicity, all stars 
with masses ^ 4OM are thought to collapse directly 
into a black hole (Heger et al. 2003) without exploding 
as supernovae. Second, the mass loss due to winds is 
much more reduced in metal-poor stars, which greatly 



helps in increasing the mass of the final remnant. 

In this paper we investigate the formation of BHs, rem- 
nants of VMS formed via stellar collisions in the very first 
stellar clusters at early cosmic times. We derive the prop- 
erties of the ensuing BH population, that can represent 
the "seeds" of the supermassive BHs that we observe in 
today's galaxies. In the following we summarize the main 
features of our model: 

• we focus on halos with virial temperatures T V i r > 
10 4 K at z > 5 after the first episode of star forma- 
tion, hence with a low , but non-zero, metallicity 

• low metallicity gas can fragment and form low- 
mass stars only if the gas density is above a certain 
threshold, n crit .z 

• halos possess angular momentum, acquired 
through interaction with their neighbors, hence 
gas collapse ultimately leads to the formation of a 
disk 

• if the forming disk is Toomrc-unstable, instabilities 
lead to mass infall instead of fragmentation into 
bound clumps and global star formation in the en- 
tire disk 

• the gas inflow increases the central density, and 
within a certain, compact, region n > n critt z 

• star formation ensues and a dense star cluster is 
formed 

• if the star cluster goes into core collapse in < 3 Myr 
runaway collisions of star form a supermassive star, 
leading to a massive BH remnant. 

We describe the physical mechanism for cluster forma- 
tion in section 2, where we discuss the conditions lead- 
ing to very compact stellar clusters in low metallicity 
(Z ~ 10~ Zq) protogalaxies. We discuss dynamical 
instabilities (Section 2.1) driving gas towards the cen- 
tre, and fragmentation occurring at high densities (Sec- 
tions 2.2 and 2.3). The properties of the star clusters, 
VMSs and BHs are described in Sections 2.4 and 2.5. 
We present our results in Section 3 and discuss the im- 
plications of our model in Section 4. 

2. BH SEED FORMATION: MODEL SET-UP 

We start by considering a halo of mass Mh, with virial 
circular velocity Vh , virial radius Ru , angular momentum 
J. The rotational support can be quantified in terms 

of the spin parameter A = Jl^l 1 / 2 / 'GM^ 2 , where E is 
the total energy of the halo. We further assume that 
a fraction ~ 0.05 of the halo mass is in a gaseous 
baryonic component than can cool and condense. The 
collapsing mass will then be vrid Mh ■ The angular mo- 
mentum of the collapsing baryons will be a fraction jd of 
the halo angular momentum, with jd/trid = 1 if the spe- 
cific angular momentum of baryons is conserved during 
collapse. In this work we focus on relatively massive ha- 
los (T v i r > 10 4 K, ensuring that gas can cool by atomic 
hydrogen) after the birth of the first stars in the Uni- 
verse. 10 4 K halos are likely build up by smaller halos 
that have already experienced a first episode of star for- 
mation. We assume that the first PopIII stars have been 
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able to affect the gas both radiatively (thus precluding 
subsequent H 2 cooling in the halo) and chemically. The 
first condition implies that, at least at the low density 
at which the collapse starts, H 2 is dissociated and does 
not contribute to the cooling of the gas (see also Figures 
9 and 10 in O'Shea & Norman 2008). The second con- 
dition ensures that metals (and dust) in small quantities 
can cool down the gas efficiently only as the gas density 
reaches a critical threshold (see Section 2.2, and Figures 
2, 4 and 5 in Smith et al. 2008). Jointly, the two con- 
ditions ensure that at first cooling is driven by atomic 
hydrogen only. 

2.1. Disk formation and mass inflow 

We describe the early evolution of the cooling and col- 
lapsing baryons in a simple way, assuming adiabatic re- 
sponse of the halo to gas cooling and disk formation (e.g., 
Mo, Mao & White 1998, Oh & Haiman 2002, Lodato & 
Natarajan 2006). As we are assuming negligible H 2 , the 
tenuous gas cools down by atomic hydrogen only until it 
reaches T gas <~ 4000 K. At this point the cooling function 
of the atomic hydrogen drops by a few orders of magni- 
tude, and contraction proceeds nearly adiabatically. 

Given the presence of angular momentum, the con- 
traction in the equatorial plane stops as the gas becomes 
rotationally supported, and the collapse ultimately leads 
to the formation of a disk. As pointed out by Mo, Mao 
& White (1998), the properties of the disk can be deter- 
mined by assuming a profile for the surface density S. 
As the mass assembles in the protodisk the surface den- 
sity increases. If for the final configuration the Toomrc 
parameter: 
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where c s is the sound speed, k is the epicyclic frequency 
and Q c is a critical value, then the disk becomes unsta- 
ble, and it develops bar-like structures. We argue that 
if the destabilisation of the system is not too violent, 
instabilities lead to mass infall instead of fragmentation 
into bound clumps and global star formation in the en- 
tire disk (Shloshman et al. 1990, Lodato & Natarajan 

2006). This is the case if the inflow rate is below a crit- 
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ical threshold M max = 2a c ^j the disk is able to sustain 
(where a c ~ 0.12 describes the viscosity) and molecular 
and metal cooling are not important. We substantiate 
further our argument in Section 2.2. In the following, we 
discuss what happens to the disk if fragmentation is (at 
least initially) suppressed. 

The classical analysis yields Q c ~ 1 for axi-symmetric 
disturbances. However a disk can develop large scale 
non-axisymmetric instabilities even for Q > 1 so that a 
limit value Q c <~ 2 is probably a better guess (Lodato 
2008). We will assume in the following Q c = 2. The fre- 
quency of disk instabilities as a function of Q c is discussed 
in Lodato & Natarajan (2007) and in Volontcri, Lodato 
& Natarajan (2008), who showed that any Q c < 1.5 un- 
derestimates the BH population, even assuming a 100% 
of BH formation in unstable disks. 

The torques induced by self-gravity can cause a rela- 
tively fast redistribution of gas within the disk: the gas 
shocks and loses angular momentum thus sinking to the 
centre of the system. As gas flows into the central region, 
the surface density in the outer part of the disk decreases 



until the Toomre parameter nears the critical value once 
again. At this point the infall stops. In this sense the 
instabilities described here are self-regulating and end up 
in a condition of marginal stability. Ultimately, the mass 
routed in the central part corresponds to the amount nec- 
essary for the outer disk to be marginally stable. Let us 
define this central mass condensation as a fraction m a 
of the halo mass, that corresponds to a baryonic mass 
m a M h . 

A marginally stable, isothermal disk has the surface 
density profile consistent with a Mestel disk, with S = 

^° (^0) ' ^ s our initial condition for the outer re- 
gion of the disk. As the mass m a M h is added to the mass 
already in place in the centre, the inner surface density 

profile steepens, S cx (j^J 1 with 7 > 1 . Mincshigc & 

Umemura (1997) studied the collapse of an isothermal, 
rotationally supported, self-gravitating collapsing disk. 
The collapsing disk develops a different structure in the 
inner and outer part, and the collapse proceeds with time 
as follows: 



S(i?) ~ < 



Sm i (£) R<R 



2/3 



R > R* 



(2) 



where the transition radius Rt r between the two profiles 
evolves linearly with time. As the transport associated 
with the development of gravitational instabilities can be 
well described in terms of an effective viscosity (Lodato 
& Natarajan 2006), we take the result of Mineshige & 
Umemura as a reference for our system and we assume 
that the inner profile steepens with 7 = 5/3 unless oth- 
erwise noted. 

The final configuration of the structure would then be 
characterized by an outer surface density profile of a Mes- 
tel disk and a central denser region: 



S(i?) 



R < Rt 



R > Rt 



(3) 



where T, in can be written as S (i?o/-Rtr ) 1 ~ 7 by imposing 
continuity at R tr . At this point the parameters of the 
final disk are So, Ro, Rtr and the fractional mass of 
the halo m a that participate to the infall. They can 
be determined by imposing the conditions of mass and 
angular momentum conservation, Q = Q c (for R > R tr ) 
and by imposing that a mass m a Mh is added to that 
already in place inside Rt r - For simplicity we assume 
that the dark matter halo hosting the protogalaxy follows 
an isothermal profile. The resulting disk properties are 
as follows: 

10m d (m d /j d ) 2 H z V h {l-m a /m d ) 3 
^0 = ) ( 4 J 



Rq — 



16ttGA 2 
2V2( Jd /m d )XR h 
(1 - m a /m d ) 



m d 



1 - 



8 ^ 3d ( Tg as 



m d Q c m d \ T, 



1/2 



(5) 
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m a M h 
4irZ R ' 



(7) 



where H z is the Hubble constant at redshift z. The 
vertical structure of the disk is determined by solving 
the equation for hydrostatic equilibrium. For a disk 
that is isothermal and self-gravitating, this implies a 
z-dependency for the density oc cosh~ 2 (h/H) where H 
is the vertical scale-height and it is equal to H(R) — 
Cg/%G'S(R). We can therefore express the inner disk 
density as: 



-5/3 

n(R,h) = no I ) cosh~ 2 (h/H(R)), (8) 



R 

,Ro, 

where no is a function of the disk parameters, So, Rq 
and 7 (see Appendix A). 

2.2. Fragmentation of the disk 

A necessary condition for the inflow process described 
in the previous Section is that star formation does not 
take place in the entire disk. If this happens, the gas that 
would otherwise flow into the central region, is consumed 
as it is converted into stars. 

Rice et al. (2005) proposed that fragmentation in a 
thin disk sets in when the gravitationally induced stress 
exceeds a critical value. Describing angular momentum 
transport in terms of the a prescription for viscous dissi- 
pation (i.e. that the torque strength can be expressed in 
terms of a, Shakura & Sunyaev 1973), then the critical 
threshold for fragmentation a c determines how much an- 
gular momentum can be transported in a steady state. 
In this sense the fragmentation boundary is ultimately 
due to the inability of the disk to redistribute angular 
momentum on a sufficiently short timescale. In the clas- 
sical analysis (Gammie 1996) fragmentation develops as 
the condition t coo i = td yn is reached. In the framework of 
Rice et al. (2005) we can expect that the same effect de- 
velops if the mass inflow from the halo induces too strong 
a stress. Lodato & Natarajan (2006) apply this frame- 
work to disks of primordial composition. By requiring 

v 3 

the mass- accretion rate from the halo Mh — rnd-fj- to be 

_3 

the disk is 
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less than the maximum value M m 

able to sustain, they argue that to avoid fragmentation 
it must be 



T 



Tgas 



< 



4a r 
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where a c ~ 0.12 (Clarke et al. 2007) is the critical value 
for fragmentation. 

We assume T gas — 4000 K, which corresponds to 
M max ~ 1O _2 M0 yr _1 . The joint conditions, m a > 
(Equation [5]) and Equation [5] then impose an upper limit 
to the virial temperature T„ r ^ 1.8 x 10 4 K. We therefore 
consider halos with 10 4 T mr iS 1.8 x 10 4 K. 

By considering these gas and virial temperatures we 
can ensure that the fragmentation threshold is not at- 
tained at least at R > R tr - However, T gas — 4000 K 
requires the absence of coolants other than atomic hy- 
drogen. If i?2 contributes efficiently as a coolant, the gas 



temperature drops, and Equation [9] implies that massive 
halos are subject to strong fragmentation. For instance, 
if gas cools to Tg as ~ 200 K, only halos with T„ r ^ 1000 
K satisfy Equation [9] At redshift ~ 20 this corresponds 
to halo masses around a few 10 5 M Q . 

The issue of H2 cooling suppression in the presence of a 
strong UV field has been addressed in a number of recent 
studies on the first galaxies. Even if a strong UV back- 
ground is not already in place, a "local" UV field can be 
established by the presence of PopIIIs in surrounding ha- 
los. Eventually, H2 can be suppressed by previous PopIII 
star formation in the smaller systems assembled to form 
the final T mr ~ 10 4 K halo (see Ciardi (2008) or Ciardi 
& Ferrara (2005) for a discussion on this issue). O'Shca 
& Norman (2008) have shown that halos with T„ r ~ 10 4 
K embedded in a UV background with strength higher 
than ~ 3 x 10~ 23 erg s _1 cm' 2 Hz" 1 sr" 1 are not sub- 
ject to efficient cooling and fragmentation in the disk (see 
their Fig. 9b) and it is only in the central region (r < 10 
pc) that Tg as drops down to ~ 10 3 K, leading to the for- 
mation of a central single massive star born in isolation. 
Halos subject to a weaker UV field show the same radial 
behavior but shifted at lower temperature. If H2 cooling 
in fact does not act on large scales, it does not affect the 
inflow process described in Section 2.1 so that material 
can still be transported to small distances by dynamical 
instabilities. 

In metal-free conditions, the gaseous density in the cen- 
tral region increases until H2 is activated (O'Shea & Nor- 
man 2008, Omukai et al. 2008, and references therein) 
and the formation of a Pop III star can proceed. On 
the other hand, if the gas has been enriched to a certain 
level, fragmentation can take place and an entire stellar 
cluster is formed (Clark et al. 2008). It is this last situ- 
ation that we want to examine in more details, i.e. the 
possibility that a cluster of stars is formed instead of a 
single massive star. 

2.3. Critical metallicity for fragmentation 

Various authors have suggested that the presence of 
a certain amount of metals is the key ingredient in or- 
der to produce efficient fragmentation (Schneider et al. 
2006; Bromm et al. 1999, 2002; Omukai et al. 2008; 
Clark et al. 2008 and references therein). Let us define 
a critical metallicity Z cr u where transition from PopIII 
to "normal" star formation occurs. Z cr n depends on the 
coolants at work. Bromm et al. (2001) simulated the 
collapse of a halo of 2 x 10 6 Af Q at different values of 
Z . They estimate Z cr n ~ 10~ 35 and show that when 
Z > Z cr n a rotationally supported disk that fragments 
vigorously can form. More recently, Clark et al. (2008) 
simulate the collapse of a rotating cloud polluted by dust. 
They show that a tightly packed cluster of protostars is 
formed in the centre. 

To estimate the value of Z cr n at which gas start to 
fragment one can compare the cooling rate A coo ; to the 
adiabatic heating rate T a d- Fragmentation requires the 
A CO oi = r ac ;. Once this happens, the temperature starts 
to decrease as the density increases until cooling stops to 
be efficient. The condition T a d < A coo / in fact constraints 
the metallicity at a given temperature and density to 
be greater than a critical value Z cr n- For an initially 
isothermal gas this condition provides a relationship be- 
tween metallicity and density. Equivalently, if the gas is 



Formation of massive black hole seeds 



5 



characterized by a fixed metallicity Z, then only those re- 
gions with density greater than a given threshold n cr it,z 
are able to cool down efficiently and fragment into stars. 
This is the condition that we apply to collapsing disks to 
determine if they can develop into stellar clusters. 

Santoro et al. (2006) studied the condition for frag- 
mentation in a low metallicity Universe assuming that 
PopIII stars are the main sources of pollution. The 
Z C rit ~ n crit.z relation depends on the ratio between the 
different species of coolants (see Figure 9 in Santoro et al. 
2006). In all cases, the critical metallicity increases with 
decreasing density untill n reaches the critical value for 
collisional de-excitation of the dominant coolants. We 
adopt as a reference the curve corresponding to PopIII 
stars in the intermediate mass range (mass range 185- 
205M Q ; however the solar abundance ratio produces a 
similar pattern as well, cfr. solid and dotted curves in 
Figure 10 in Santoro & Shull 2006). We discuss how our 
results depend on the specific Z cr n — n cr it,z relation in 
section 3.2.2. 

The metallicity of gas in a given halo most probably de- 
pends on its mass and on the redshift. The cosmic metal- 
licity history (MEH) has been investigated both from 
an observational and a theoretical prospective (Scanna- 
pieco et al. 2003, Tornatore et al. 2006, Savaglio 2006, 
Savaglio et al. 2005, Prochaska et al. 2003, Prochaska 
et al. 2007, Kulkarni et al. 2005, Li 2007). The ob- 
served mean metallicity decreases with redshift, but the 
rate of decrease depends on the type of sources studied. 
From studies of QSOs-DLAs, Z/Zq scales as 10~ f3z with 
(3 = 0.36 (Li 2007), while measurement based on GRBs 
point towards a shallower decline, so that at z ~ 3 — 4 
already Z = lO" 1 ^ (Savaglio 2006). The differences 
in the slope and in the normalization are often ascribed 
to different histories of metal enrichment in QSOs and 
GRBs hosts (Li 2007 and references therein). 

We model the MEH based on fits presented in Li 2007, 
that we refer to for additional details. In our reference 
model we adopt (3 = 0.36, that is Z/Z Q = 0.35 x lO" ' 362 . 
The case (3 = 0.18, leading to Z/Z Q = 0.63 x lO" 0182 , 
is discussed in section 3.2.3. We further allow a logarith- 
mically uniform scatter in Z of Alog(Z) = 1.5, based 
on the observed scatter in Z / Zq from measurements of 
the QSO-DLAs. We wish to stress that our treatment 
of MEH is highly simplified. In a forthcoming paper we 
will determine self-consistently the evolution of the MEH 
from stellar winds pollution. 

2.4. Stellar cluster formation 

The central star cluster forms within the region where 
gas at a given metallicity has reached the critical den- 
sity for fragmentation. Given a specific density profile 
(Equation [2]) the requirement n > n cr i t ,z translates into 
a condition R < Rsf, where we define Rsf as the radius 
where the density in the plane h = 0, reaches the value 
Merit, z for the gas metallicity Z, that is: 

n(R = R S f, h = 0) = n crittZ - (10) 

As long as n(Rt r ,0) < n cr it,z, n(R,h) is defined by 
Equation [8] and the radius within which gas fragments 
and stars form can be written as: 

l/ 7 

E Ro_ / ttG 

R tr V 2nm H n criti z 



Rsf = Rtr 



(11) 



where the metallicity dependence is included implicitly 
in n C rit,Z (see Appendix A), /i is the molecular weight 
and trig is the proton mass. The half mass radius can 
be expressed as a function of the cluster radius in a very 
simple way, as 



(12) 



For a given 7 the mass in stars of the cluster can be 
calculated as 



M cl = e SF 2n 



E(R)RdR 



: £sf 
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/ ttG 
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(13) 



where esf is the fraction of gas converted into stars. We 
assume e$F — 0.25, consistently with the star forma- 
tion efficiency in the low redshift Universe (Lada & Lada 
2003). Given the metallicity, Z, it is now possible to de- 
termine the extent of the region inside which star forma- 
tion is allowed. The cluster properties for a given halo, 
are therefore uniquely described for any given metallicity. 

At high metallicities, the case n(R tr ,0) > n cr i t: z be- 
comes common. This implies Rsf > Rtr, that is, star 
formation takes place in the region of the disk where 
no inflow is taking place. In other words, once the gas 
reaches Rsf, it is no longer able to be routed efficiently 
in the inner region as star formation starts to consume 
gas. In the remainder of the paper we will assume conser- 
vatively that if Rsf > Rtr no cluster formation occurs, 
as stars are formed in the disk rather than in the central 
compact region. 

2.5. Runaway instability of the central cluster 

The conditions under which "mass segregation insta- 
bility" can occur have been investigated in a series of 
papers, focusing on clusters in the present-day Universe. 
A successful core collapse requires that the core collapse 
time be less than the main sequence lifetime of the most 
massive stars (mass losses from supernovae expand the 
core and increase interaction times) . The main sequence 
lifetime of massive stars asymptotes to about 2.5 Myr 
because all stars go off the main sequence when they 
have consumed about 15% of their hydrogen, and for 
high-mass stars with luminosities approaching Edding- 
ton, L cx M (not M 3 ' 5 as is the case for lower masses) 
and then the lifetime cx 0.15M/L ~ const. 

With a Monte Carlo code, Giirkan et al. (2004) found 
that typically the mass of the collapsing core is 10~ 3 
times that of the entire cluster. Similarly, Portegies 
Zwart and McMillan (2002) related the mass of the VMS 
with the parameters of the cluster taking into account 
both numerical simulations and analytical arguments. 
Portegies Zwart & McMillan (2002) have shown that core 
collapse occurs on a timescale 
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where In Ac is the Coulomb logarithm of dynamical 
friction (Binney & Tremaine 1987). The same simula- 
tions also find that the mass of the VMS can reach val- 
ues as high as 10 3 M Q (see also Freitag et al. 2006b). 
The final mass of the VMS however depends on complex 
phenomena related to both the dynamics and hydrody- 
namics of the collisions. At solar metallicity the growth 
of a star with mass greater than ~ 100M© can be highly 
problematic as mass loss occurs both during the main 
sequence phase and at the end of the evolution when the 
star collapses into a BH. The growth of a VMS should 
in principle be much more efficient for metal-poor stars 
as in this case as mass loss should be strongly reduced. 
Recent models of stellar evolution at Z ~ 10 -5 Z§ have 
shown that mass loss due to stellar winds during the main 
sequence phase is almost unimportant and that the main 
contribution to the reduction of the stellar mass is due 
to the effect of rotation. This can reduce the mass of 
the star by a factor of order 2-4 (Meynet et al. 2008). 
If the final mass achieved by the VMS is greater than 
~ 260M Q , then after the main sequence it collapses into 
a BH retaining most of its mass (Heger et al. 2003). 

To estimate of the final mass of the VMS we follow 
the treatment outlined in Portegies Zwart (2002). We 
assume that the mass of the BH seed, Mbh, corresponds 
to the final mass of the VMS: 

M BH = m» + 4x lO" 3 M cl f c In \ c In ' ' ' ' " ^ 



V tec J 



(15) 



Here tMS = 3 Myr is the main sequence lifetimes 
of massive (> 40 M®, Hirschi 2007) stars, to* is the 
initial mass of the seed star that experiences runaway 
growth and f c is a factor used to calibrate the ana- 
lytical expectation with direct numerical simulations. 
Portegies Zwart & McMillan (2002) find f c = 0.2 for 
ln(Ac) = ln(0.1M c //(m)). We adopt the same values as 
Portegies Zwart & McMillan (2002) for both parameters, 
further assuming to* = (to) = 10M Q consistent with the 
characteristic stellar mass in 10 7_8 Mq halos at redshift 
~ 10 (see Figures 2 and 3 in Clarke & Bromm 2003). 

It is important to stress that Equation [15] provides an 
upper limit to Mbh- However, once the seed is born, 
the BH is still embedded into a dense cluster of stars. 
Even if its growth has been limited before the collapse of 
the VMS, the remnant can still gain mass by accretion 
of stars (for the growth rate of a massive BH hosted in a 
cluster of stars see for example the models of Marchant 
& Shapiro 1979). If the combination of star formation 
and stellar feedback of the cluster stars do not deplete 
of gas the inner region, this surviving gas can supply an 
additional reservoir of material for BH growth. Both this 
processes can contribute to bring the BH mass to values 
as high as Equation [T51 would suggest. 

Clusters forming in unstable disks might have some 
degree of rotation. In this case the gravogyro instability 
(Inagaki & Hachisu 1978) could contribute to accelerate 
the core collapse. The gravogyro instability is believed to 
occur in systems that exhibit a radial gradient of the an- 
gular speed. In analogy with viscous transport, angular 
momentum is transfered outwards. The core of the star 
cluster then contracts because of a deficit in the centrifu- 
gal force. Ernst et al. (2007) show that the gravogyro 
instability occurs in clusters with equal-mass stars, but 
in systems with two-mass components the effect of rota- 
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Fig. 1. — Cluster properties as a function of the spin parameter 
A = J\E\ 1 / 2 /GMt (i.e., the fraction of halo support given by 
rotation) for 2 critical densities for fragmentation, n cr a = 10 3 
cm -3 (solid line), and n cr n = 10 4 cm -3 (dashed line). Curves are 
truncated when Rgp equals Rtr- 
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Fig. 2. — Cluster properties as a function of metallicity for 3 
spin parameters: A = 0.018 (solid line), A = 0.01 (dashed line), A = 
0.015 (dotted line). Each curve begins at low metallicity when halos 
starts star formation and ends when Rsf equals Rtr', this happens 
at different critical mctallicitics for different spin parameters. Here 
the Z cr i t — n cr it t z relationship follows the intermediate PopIII 
mass case in Santoro et al. 2006. 

tion seems negligible, as mass segregation and rotation 
compete in leading the evolution of the stellar cluster. 
For this reason we neglect its contribution in this work. 

2.6. Cluster Properties 
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Various cluster properties: size (Equation 11), half- 
mass radius (Equation 12), mass (Equation 13) are 
shown in Figure [I] as a function of the spin parameter 
for two representative values of the critical density for 
fragmentation. We also present BH masses (Equation 
15), where we have further imposed the condition t cc < 3 
Myr. 

As A increases, less mass inflows within R tr , and Rsf is 
reached at smaller radii at a given n crit . Therefore Rsf 
and M c i decrease with increasing A, contrary to the disk 
size i?o- Note how at low critical densities (correspond- 
ing to high metallicities) only a very small fraction of sys- 
tems can undergo core collapse, as on the one hand clus- 
ter formation is suppressed for large A (as Rsf > Rtr), 
on the other hand clusters are too massive and large at 
small A for fulfilling the condition t cc < 3Myr. Clusters 
undergoing core collapse have preferentially large spin 
parameters, within the region where cluster form. Fig- 
ure [2] shows cluster sizes, masses, core collapse timescales 
(Equation 13) and BH masses as a function of metallicity. 
As the metallicity increases, the critical density for frag- 
mentation decreases and cluster formation is eventually 
suppressed when n(Rt r ,0) > n crvt .z (Rsf > Rtr)- At 
metallicities below 10 -5 Zq fragmentation is impossible 
for the specific choice of the Z cr n — n C rit,Z relationship 
(intermediate PopIII mass case in Santoro et al. 2006). 

The parameter space (virial temperature, spin parame- 
ter) where the multiple instabilities are efficient is shown 
in Figure [3] Here we select halos with T„j r > 10 4 K at 
z = 12, and derive the disk and cluster properties, as- 
suming a single critical density for fragmentation (from 
10 3 cm~ 3 to 10 4 cm~ 3 ). The higher the critical density for 
fragmentation (i.e., the lower the metallicity, see Section 
2.3) the more compact are the clusters, and the shorter 
is t cc . When n crity z < 10 3 cm -3 no clusters can undergo 
core collapse in less than 3 Myr. When n cr i t ,z > 10 4 
cm~ 3 all forming clusters undergo core collapse in less 
than 3 Myr. 

3. RESULTS 

3.1. Summary of the general procedure 

We first summarize the procedure taken in order to 
determine the properties of the BH seed population. 
We calculate the mass of halos that at redshift z cor- 
respond to virial temperatures 10 4 T v i r 5^ 1.8 x 10 4 
K (Barkana & Loeb 2001) and we determine their fre- 
quency using a modified version of the Press & Schechter 
formalism (Sheth & Tormen 1999) in a WMAP5 cosmol- 
ogy (Dunkley et al. 2008). To each halo we assign a value 
of the spin parameter, A, extracted from the probability 
distribution found in the Millennium simulations (Bett 
et al. 2007): 



/>Mo;;A) A (A 



exp 



3/C 



(16) 



where Ao = 0.0043 is the peak location, ( — 2.509 and the 
normalization reads A = 3 In 10C'~ 1 /r(C), with T being 
the gamma function. This set of assumptions allows us 
to calculate the initial disk properties, So, Ro, Q c , Rtr 
and m a . 

We assign to each halo a metallicity, Z by extrapolat- 
ing at higher redshift the fit to the observational con- 
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Fig. 3. — Parameter space (virial temperature, spin parameter) 
for cluster and BH formation. Here we select halos with T v i r > 10 4 
K at z = 12, and derive the disk and cluster properties, assuming 
a single critical density for fragmentation (top-left: 10 4 cm -3 , top- 
right: 2 X 10 4 cm -3 , bottom-left: 2.5 X 10 4 cm -3 , bottom-right: 
3 X 10 4 cm — 3 ; for the chosen Z cr u — n cr i t .z relationship these den- 
sities correspond respectively to log(Z/ZQ)=— 4.2;-4.3;— 4.4;-4.5). 
The black shaded area shows the range of temperatures and spin 
parameters where disks are Toomre unstable and the joint con- 
ditions, m a > (Equation [SJ and Equation [9] are fulfilled. The 
lighter shaded area selects the systems where Rsf < Rtr- The 
hatched area picks the subsample of clusters where t cc < 3Myr, 
where VMSs and BH seeds can form. 



straints of the MEH (Z oc 10~^ z ), taking also into ac- 
count the observed metallicity scatter. We then calculate 
ncrit, z from a given Z cr a — n cr i tt z relation. If a proto- 
galaxy has Q < Q c we determine the properties of the 
stellar cluster (Rsf, M c i and Rh)- We then check if the 
cluster can develop runaway instability via Equation ll41 
and we select the systems where t cc < 3Myr. For these 
unstable clusters we calculate the expected mass of the 
seed BH from Equation [15j In Table [T] we summarize all 
the different cases we describe in the next sections. 



Run 


MEH 




Alog(Z) 


A 


0.36 


2 


1.5 


B 


0.18 


2 


1.5 


C 


0.36 


2 





D 


0.36 


1 


1.5 


E 


0.36 


3 


1.5 



TABLE 1 

List of the simulations: label of the run, power-law index 
of the metallicity dependence by redshift, popiii yield 
curves (1,2 and 3 refer to lower, central and upper 
curves of Figure 10 in Santoro & Shull (2006) 
respectively), allowed variation (logarithmically 
uniformly distributed) in z. 



We now discuss the properties that nuclear clusters 
possess at birth and the resulting BH seed population. 
We start describing our reference model, A, in Section 3.1 
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Fig. 4. — Mean cluster masses (upper left panel), radii (upper 
right), core collapse timescales (lower left) and BH masses (lower 
right) as a function of rcdshift for model A. The orizontal line in 
the lower left panel marks the critical core collapse timescale for 
VMS formation. The dashed curves mark the 1— cr scatter. 
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Fig. 5. — Top: Cluster mass function for model A integrated over 
all redshift. The shaded area corresponds to the subpopulation of 
clusters able to form BH seeds. Bottom: BH seed mass function 
integrated over all redshifts. 



while in Section 3.2 we discuss how our results depend 
on PopIII stars metallicity patterns and on the rate of 
metal enrichment of the Universe. 

3.2. Model A 

The first halos reach the critical metallicity for frag- 
mentation at rcdshift ~ 14. This is also when the first 
stellar clusters form. The first clusters have masses of 
the order of 10 5 M Q and Rh ~ 0.5 — 1 pc; in such com- 
pact clusters, core collapse starts early ((t cc ) ~ 0.1 Myr 
at z=14). At later cosmic times, the average gas metal- 
licity increases, so that the critical density for fragmenta- 
tion decreases. A lower critical density implies that R$f 
increases (Equation ITT]) . lengthening the core collapse 
timescale. Therefore, clusters form less concentrated and 
more massive and their core collapse timescale continues 
to increase with decreasing redshift. This behavior is 
evident in Fig. [4] for the mean quantities ((M c i), {Rh}, 
(t cc ) and (Mbh))- The mean seed mass as a function of 
redshift is shown in Figure [4] (lower right panel) . Unlike 
(M c i) and (Rh), (Mbh) shows no redshift dependence as 
the increase in M c i and t cc compensate (see Equation II 5[) 
leading to a roughly constant Mbh- 

In Figure [5] the entire cluster mass function (solid his- 
togram) is compared to the mass function of systems 
able to form a BH seed (dashed histogram). The clus- 
ters that do not form BH seeds are the very last to form, 
when the metallicity of the Universe is already signifi- 
cant. The BH mass function is also shown in Figure 
The mass function is peaked at ~ 1000M Q with a long 
tail at low masses, and a steep drop at high masses. 

Fig. [8] shows the fraction, /bh, of halos hosting a BH 
seed. Seeds start to form at z ~ 14 in coincidence with 
the first, very compact, stellar clusters. The typical core 



collapse timescale increases with cosmic time, (t cc ) ~ 3 
Myr at z ~ 9, and by z ~ 8 /bh drops rapidly. From this 
point on even the most concentrated and least enriched 
disks are unable to create central stellar concentrations 
with t cc < 3 Myr and seed formation is completly sup- 
pressed. Figure [5] also shows the integrated comoving 
mass density of seeds p see d- At z ~ 10 the mass den- 
sity saturates at a value of ~ 300 M Q Mpc -3 . This mass 
density should be considered a lower limit for the total 
black holes mass density /?bh as we are completely ne- 
glecting BH growth after seeds formation. We will dis- 
cuss the implications of this formation mechanism on the 
evolution of the supermassive black hole population in a 
forthcoming paper. We note however that this seed den- 
sity is similar to that expected from PopIII star remnants 
(roughly a factor of 3 larger, Volonteri, Haardt & Madau 
2003), and we therefore expect that most observational 
constraints can be fulfilled at the same level. 

3.3. Impact of the uncertainties on metal enrichment 
onto the seed population 

We now discuss how our results depend on our choice 
of parameters. We compute the seed population for BHs 
formed over cosmic time varying (i) the history of metal 
enrichment (models B and C) and (ii) the ratio of the 
coolants that determine the Z cr n — n cr n.z relationship 
(models D and E). 

Mass functions of protogalactic nuclear clusters for all 
models defined in Table 1 are shown in Figure 12; Cluster 
stellar masses, M c i, span 2 orders of magnitude, between 
10 4 - 10 6 M o , with a peak around ~ 10 5 M Q . Figure 
[2] also shows the distribution of core collapse timescale 
of all clusters. The dashed vertical line at 3 Myr marks 
the limit for VMS formation. As discussed in the previ- 



Formation of massive black hole seeds 



9 




Fig. 6. — Top: Cluster mass functions integrated over all redshift. 
Bottom: distribution of core collapse timescales, integrated over 
all redshifts. Lines color and style as follow. Solid: black = A. 
Dotted: thin =B, thick = C. Dashed: thick = D, thin = E. Model 
A is shown in all panels as a reference. 

ous section, clusters satisfying the condition for the on- 
set of runaway instability are clustered at small masses 
and radii: typical masses of runaway unstable clusters 
are around a few 10 5 Mq and typical radii are ~ 1 pc 
(see Equation [14|) . This naturally points towards the 
very first clusters: the first systems that form are indeed 
those that more easily can give birth to BH seeds. Con- 
sequently, these are also the most metal-poor clusters, so 
that our picture is consistent with requiring that VMS 
can more easily grow in low metallicity environments. 

As times goes on, both (M c i) and (Rh) grow: as a 
consequence (t cc ) increases. Even if clusters continue to 
form, BH seeds cannot be created any longer (cfr. Fig- 
ure [7]). The mass functions of BH seeds are shown in 
Figure [3 The seed mass distributions show a character- 
istic shape with a peak at a few 10 3 M© and a long tail at 
lower masses with very little redshift dependence. This 
general picture is valid for all models; we now discuss in 
turn the dependencies on specific model parameters. 

3.3.1. Changing MEH 

The MEH is one of the most uncertain parameters. We 
explore two extreme cases in models B (/3 = 0.18, blue 
dotted curves in Figure 6 and onwards) and C (f3 = 0.36 
and Alog(Z) = 0, green dotted curves in Figure 6 and 
onwards). The MEH determines, together with the cho- 
sen Z cr it — n cr it,z relation, when seeds form. The dura- 
tion of the BH seeds formation epoch is indeed given by a 
combination of the assumed metallicity spread, A log(Z), 
and of the slope of the Z(z). Once Z increases over 
Zcrit.min, fragmentation is activated in a more extended 
region of the disk, the inflow is reduced and runaway in- 
stability can not proceed efficiently. The seed formation 
epoch is therefore longer either if Z has larger A \og(Z) 



Fig. 7. — Top: Mass function of seed BHs, integrated over all 
redshifts. Bottom: mean BH mass as a function of redshift. Lines 
color and style as in Figure [6] 



or if metal enrichment is rather inefficient. 

Model B has the most efficient metal enrichment, and 
BHs appear already at redshift 30. On the other hand, 
BH formation is also suppressed very early, at z ~ 18. 
At this early cosmic epoch very few halos were massive 
enough for efficient atomic line cooling, thus leading to 
a comoving seed mass density of only a few M . As 
previously noted for model A p see d does not necessary 
coincide with /jbh as we are neglecting seeds growth. In 
model B seeds form at higher redshift, and mass can be 
built up for a longer period, likely increasing BH masses 
by accretion. 

Model C has the same redshift dependence of the MEH 
as our reference model A, but we assume no scatter. A 
null A log(Z) (model I) produces a short burst of seeds 
very concentrated in time, as seed formation is allowed 
only for a very sharp range of Z. Clusters and BHs form 
m a burst at z = 11 - 12 when Z ~ 10~ 5 (the minimum 
in the Z cr i t — n cr i ti z relation). This burst is very efficient, 
with a high fraction of halos undergoing cluster and BH 
formation, and the resulting p see d ~ 300M Q Mpc~ 3 is 
similar to our reference case. 

3.3.2. Changing Z cr it — n cr it,z 

We explore the effect of different yield patterns for 
PopIII stars on our results in models D (blue long-dashed 
curves in Figure 6 and onwards) and E (green long- 
dashed curves). The adopted Z cr n — n C rit,z relation de- 
fines together with the MEH, when seeds are born. The 
epoch at which BH seeds start to form is related to the 
redshift at which the first systems reach the minimum 
possible critical metallicity Z crit min ~ 10~ 5 Z & (see the 
minimum of the curves in Fig. 10 of Santoro & Shull 
2007). 

If fragmentation is allowed at lower metallicities 
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Fig. 8. — Top: fraction of halos hosting a BH seed as a function 
of redshift for models in Table 1. Bottom: Co-moving mass density 
of the seed BH for the same models. Colors and lines styles as in 
Figure [6] 



(model D) clusters begin to form earlier. With increas- 
ing cosmic time, Rsf increases to Rt r , thus precluding 
the formation of a seed. In model D this happens al- 
ready at z ~ 12, while BH formation proceeds all the 
way to z — 5 in model E. Consequently, the seed mass 
density p see d ranges from ~ 100 to 300 M Q Mpc~ 3 at 
z = 5, with model D, that forms seed early on, having 
the lower seed density. 

4. DISCUSSION 

We described a model for BH seed formation as result 
of multiple successive instabilities. On a large scale, grav- 
itational torques in Toomre unstable primordial disks 
can pile up a significant amount of gas in the central 
region of high redshift halos (e.g., Lodato & Natarajan 
2005). We focus on relatively massive halos at high red- 
shift (T V i r > 10 4 K, z > 10) after the very first stars 
in the Universe have completed their evolution. This set 
of assumptions ensures that (i) atomic hydrogen cooling 
can contribute to the gas cooling process, (ii) a UV field 
from the first stars has heated the gas, precluding H 2 
cooling except in the highest density regions , and (iii) 
the gas inside the halo has been mildly polluted by the 
first metals. The second condition implies that, at least 
at the low density at which the collapse starts, H 2 is dis- 
sociated and does not contribute to the cooling of the gas 
(see also Figures 9 and 10 in O'Shea & Norman 2008). 
The third condition ensures that metals (and dust) in 
small quantities can cool down the gas efficiently only as 
the gas density reaches a critical threshold. Jointly, the 
two conditions ensure that cooling is driven at first by 
atomic hydrogen only. 

At low metallicity (Z ~ 10 _5 Z Q ), fragmentation and 
low-mass star formation can occur only where the den- 
sity is above a metallicity dependent critical threshold 



(Z cr it — n cr it.z, Santoro & Shull 2006), corresponding to 
the central densest region of the protogalaxies. Within 
this limited region where star formation takes place effi- 
ciently, very compact stellar clusters form. The typical 
stellar masses are of order 10 5 M© and the typical half 
mass radii ~ 1 pc. Eventually, a large fraction of these 
very dense clusters undergo core collapse before stars are 
able to complete stellar evolution and a VMS can form 
(e.g., Gurkan et al. 2004, Portegies Zwart et al. 2006). 

Clusters unstable to runaway collisions are always the 
first, less massive (~ 10 5 M Q ) ones. As the metallicity 
of the Universe increases, the critical density for frag- 
mentation decreases and stars start to form in the entire 
protogalactic disk so that i) accretion of gas in the centre 
is no longer efficient and ii) the core collapse timescale 
increases. As a result less and less compact clusters form, 
and less of them are subject to rapid core collapse. 

We computed the properties of the BH population for a 
set of models, in dependence of various parameters. The 
epoch of seed formation is determined by the time at 
which gas in the centre of the halo can start to fragment. 
This redshift depends on the metal enrichment history 
and on the exact shape of the Z cr n — n cr i ty z relation. 
If metal pollution is very efficient and the Universe was 
enriched early, the BH formation epoch ends very early, 
when only a few halos were massive enough for efficient 
atomic line cooling. The mass and number densities of 
BH seeds are consequently very low. The fraction of 
halos hosting a BH seed depends also on the fraction of 
unstable disk (hence, the critical Toomre parameter, Q c ) 
and on the frequency of halos with Z ~ Z cr it, m in- De- 
creasing Q c has a twofold effect on the efficiency of BH 
seeding. First, as already noted by Lodato & Natarajan 
(2007), it decreases the number of bar- unstable disks as 
the Toomre criterion is satisfied for higher surface den- 
sities (requiring very small spin parameters). Addition- 
ally, the higher surface density for bar-instabilities im- 
plies smaller unstable regions (Rtr), but a larger Rsf at 
a fixed particle density. Cluster formation is therefore 
truncated at A < 0.004 for Q c — 1, leading to inefficient 
cluster and BH formation. The slope of the inner density 
profile, 7, affects /bh as it affects the fraction of runaway 
unstable systems that can form a compact cluster. A 
lower 7 produces shallower, more extended clusters with 
longer core collapse timescales (Equations [11] and [i~4")) . 
For reasonable choices of 7 > 4/3 we find that the effi- 
ciency of the BH formation is within a factor of 2 of our 
reference model. 

Most of our assumptions have been quite conservative, 
but still the population of seeds is comparable to the case 
of Population III star remnants discussed, for instance in 
Volonteri, Haardt & Madau (2003). The fraction of high- 
rcdshift galaxies seeded with a BH is about a factor of 
10 below the direct collapse case presented in Volonteri, 
Lodato & Natarajan (2008), where a seed was assumed 
to form with a 100% efficiency whenever a protogalaxy 
disk was Toomre unstable. An estimate of the degree 
of agreement between the evolution of the population of 
seeds that we have calculated in this paper and observa- 
tional constraints requires dedicated models, including 
the mass growth of the BHs after their formation, and 
how galaxy and BH mergers influence the population. 
We will present such models in a future paper. 
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APPENDIX 



DISK AND CLUSTER STRUCTURE 
Disk parameters 

Disk parameters Sqi Ro> Rtr arid m a have been summarized in Section 2.1. We only insert here the complete 



expression for R tr with the 7 dependency 



2 - 7 m a M h 
7 - 1 2^2,0^0 



As stated in Section 2.1, the particle density of the inner disk can be written as 



n(R,h)=n Q ^-(^-) rush -!/,///! //n (A 2) 
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1 


2 


3 


a 


0.03517305 


0.0317305 


0.0317305 


b 


-0.582132 


-0.572132 


-0.62132 


c 


-4.1 


-2.75 


-1.1 



TABLE 2 

List of the values of a, b and c for fitting different Z crit — n crit curves. 

where H can be computed imposing hydrostatic equilibrium and it is found out to be H = , Cs (see for 

example Oh & Haiman 2002). To simplify the notation we have defined f(R) = Rp/ RtrjRtr/R) 1 ■ We relate no 
and Eo by imposing that the surface density calculated starting from Equation IA2| follows the profile described by 
Equation 3: 



^of{R)^^m u n f{R) 



cosh 2 (z/h)dz 



- 2/j,m u nof(R)H = V2c 



n nm u f(R) 
irG 



Resolving for no we find: 



no 



Eg nGf(R) 
c 2 s 2fim u 



(A3) 
(A4) 



Star formation radius 

Once the density profile of the disk is defined, the star formation radius can be determined by imposing n(R = 
Rsf, h = 0) = n cr it,z- We allow clusters to form only if star formation is induced inside Rt r and not in the external 
part of the disk. The condition for this to happen is that nt r = n(R — Rt r , h — 0) < n cr %t y z- Rsf is then computed 
by imposing: 



"n 



Rq f Rtr \ 7 _ 
Rtr \RSF ) 



n C rit,Z- 



Inserting the expression for Uq 



Resolving for Rsf 



TrGSg {RoYfRtrV 
I I I — — I = n crit z 



2nm u c 2 \R tr J \ Rsf ) 



(A5) 



(A6) 



Rsf = Rt 



p 2/5 
■ ^tr 



TTG 



c s R tr y 2fim u n crlt ,z 

T,qRq / 



ttG 



i/ 7 



3/5 



2^m u n cri t,z 



(A7) 



where in last expression we have inserted 7 = 5/3 explicitly. 

As stated in Section 2.3, n crit z depends on the metallicity of the gas. The curves in Figure 10 of Santoro & Shull 
can be fitted by the expression 



log(Z/Z ) = a log 2 {n C rit, z) + b\og{n cr it,z) + c 
where the values of a, 6, and c have been calculated for the three curves and are reported in Table [2] 



(A8) 



The minimal critical density for fragmentation can be computed from Equation I A81 Inserting the expression for the 
critical density into Equation I A7I one finds 



i iz, \ 7 - 1 i /^ 2 -7 rn a M h 
log(RsF) = log 

7 

1 

1 

7-1 



log 



7 - 1 2^S i? 

•£ R \ 1, / ttG 

- - log 

2 \ 2[im u 



-log(n criti z) 



7 



[log 



2-7 
7-1 



A] 



Formation of massive black hole seeds 
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-(b- CWl + D'\og(Z/Z Q ) 
7 v 



7-1 



7 



[log 



2-7 
7-1 



A] 



1 



-[B- Cy/1 + D/3zj (A9) 

where in the last expression we have first expressed n cr i tj z as a function of Z and then we have inserted the redshift 
dependence of the metallicity 1 Z. All constants except 7 have been included in A, B, C, C, D and D'\ 



A = \og 



D- 



4a 



/ m a M h 
\2ttE Ro 

log 



■log 



/ ttG 

\ 2(j,m,i 



4a V b 2 



C = C'VT+D r 6 
4a 



D' 
D 



b 2 - 4ac 
D' 



Finally for 7 = 5/3 



\og{R SF ) =h + k 2 ^Jl + k 3 pz 



(A10) 
(All) 



where all constants are collected in k\ = ^{A — 0.3) + 3/5£>, k 2 = 3/5C and fc 3 = D. 

Black hole masses 

Black hole masses depend on the stellar mass of the clusters and on the half-mass radius, via the core collapse 
timescale. 

Cluster masses and timescales to core collapse can be estimated coupling Equations I A7I and IA1 II to Equation 12: 



/•Rsf 

Jo 

= 2nesF^oRoR] 
= 2nesF^oRo 



M r l = 



7?2-7 
37-1 ^SF 

2-7 
2-7 m a M h 



7 - 1 2ttEo.Ro 



7-1 ^2-7 



SF 



-esF 



2- 7 
3-3 



(A12) 



2-7 c s 



ttG 



1/5 



m a M h 
7-1 



Gtt 



ttG 



-m a M h 



4/5 



(A13) 



where in the last line we assume 7 = 5/3. 

The core collapse timescale is defined in Equation 1141 The half mass radius can be simply expressed as Rh = 
2 1 /(t- 1 ) R sf . We use Equation All to express the dependence of t cc on M c i as: 



tc 



t R 3 / 2 A-f 1 / 2 

T oR h M cl 

2-7 



r 2^R^Ml( 2 



t 2 2 7-1 



9-rrV R , 1 /(2-7) 



3/2 . 3 1-t 



7 



: t 



2 9/4 



1 



3/2 



3 

-m a M h 



_6TrY, R e 3 SF 

1 The last expression refers to the mean Z only, i.e. it is rigorously appropriate only to model I 



(A14) 
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where To is the normalization of Equation [14] and its value in Myr, M Q and pc is 0.019. In the last Equation we assume 
7 = 5/3. 

A/bh can now be found by inserting Equation I A14I into Equation [TS] 



where g is defined as 



M BH = m* + 4 • 10~ a f c In \M cl In 



9 



■ Mr, In 



TO 



Ml 



27re| F 7 I]oi?o 



-3/2 



m a M h 
7-1 



3 

2 2-7 



(A15) 



(A16) 



From Equation IA15I we can argue that Mbh first grows with increasing M c \ , it attains a maximum value and then 
it decreases, as more massive clusters have longer core collapse timescales with t cc oc M^. This last proportionality 
results from the dependence of Rsf (and consequently Rh) on M c /. 



